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In 1937 Bailey and Martyn proposed a theory of ionospheric cross modulation along 
the following lines. Radiation emitted by the disturbing transmitter is strongly absorbed 
in a region of the ionosphere where the wanted wave is refracted and attenuated. Absorp- 
tion of the disturbing radiation raises the temperature of the ionosphere. If the amplitude 
of the disturbing radiation varies with time, so will the resulting temperature increase and 
hence the attenuation suffered by the wanted wave. In the simplest experimental arrange- 
ment the wanted wave is initially unmodulated and the disturbing radiation is modulated 
at a definite audio frequency. In traversing the ionosphere, the wanted wave becomes 
modulated at the frequency of the disturbing radiation and also at twice this frequency. 

The Bailey-Martyn theory makes three quantitative predictions concerning the trans- 
ferred modulation. The first concerns its dependence on the power of the disturbing radi- 
ation, and involves only t he assumption that the effect is small enough to warrant a linearized 
description. The remaining two predictions concern the dependence of the transferred 
modulation and its phase lag on the modulation frequency. These predictions involve 
the additional assumption of a unique relation, valid for all values of the modulation fre- 
quency, between the mean thermal energy and the mean collision frequency of the electrons. 
A relation of this kind can exist, however, only for modulation frequencies much less than 
150 c/s. Since the experimental range of modulation frequencies extends to nearly 
1500 c/s, the foundations of the Bailey-Martyn theory need to be reconsidered. 

This paper describes a microscopic theory of ionospheric cross modulation. The 
velocity-distribution function for the electrons enters explicitly into the theory; it is deter- 
mined by a differential equation whose form depends on the modulation frequency. Knowing 
the distribution function, one can calculate the absorption coefficient. The form of the 
predicted absorption coefficient depends on the assumed form of the electron-molecule 
interaction law. By numerical methods, we have calculated the transferred modulation 
and its phase lag as functions of modulation frequency for a few simple interaction laws. 
The calculations show that the effects of departures from a maxwellian velocity distribution 
are indeed significant. The predictions are sensitive to the assumed form of the electron- 
molecule collision law. Although the present theory is still highly idealized, the results 
obtained suggest that further theoretical and experimental refinements could lead to an 
experimental determination of the electron-molecule interaction law in the J) region. 



1. Introduction 
1.1. The Bailey-Martyn Theory 

The passage of a radio wave through a region of 
the ionosphere in which it is partially absorbed 
raises the electron temperature in the region and 
hence changes its absorption coefficient. The ab- 
sorption process is accordingly nonlinear: the attenu- 
ation suffered by each of several radio waves passing 
through the same region depends on that suffered 
by all the others. This is the physical basis of 
Bailey and Martyn's [1934] theory of ionospheric 
cross modulation. 



1 The first five sections of this paper first appeared in slightly different form in 
a report dated March 14, 1959, prepared for the U.S. Air Force. A much more 
-complete discussion of the Bailey-Martyn theory and its modifications (sec. 1) 
and of the kinetic theory of a plasma in an alternating electromagnetic field 
(sections 2, 3, 4) has subsequently heen given hy Ginzburg and Gurevich (1960), 
.along with an exceptionally complete bibliography, but the main results of the 
present paper (sections 5, 6) do not appear to have been published previously. 
An approach very similar to that used in sections 2, 3 has been outlined by 
Bayet, Delcroix, and Denisse [1957], but their treatment of the distribution func- 
tion differs from that of sections 4, 5. We are indebted to the referee, Dr. Julius 
Cahn, for bringing this paper to our notice. 



In a typical cross-modulation experiment an 
initially unmodulated carrier wave (the "wanted" 
wave) traverses a region in which a second wave at a 
different carrier frequency (the "disturbing" wave) 
is heavily absorbed. The disturbing wave is ampli- 
tude-modulated at a definite frequency. On recep- 
tion, the wanted wave is found to be modulated at 
the same frequency and also at twice this frequency. 

The main predictions of the Bailey-Martyn theory 
do not depend on a detailed description of the 
physical processes involved but follow from a few 
general assumptions, the first of which is that the 
phenomenon admits a linearized description. The 
linearized energy-balance equation for the electron 
gas has the form 



de T 

tc kI - t 

Here e denotes the difference 



MT-T m ) 



(1.1) 



(1.2) 
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between the mean thermal energy of the electrons 
and the mean thermal energy of the molecules, k 
denotes the absorption coefficient of the medium, 
/ denotes the intensity of the ambient radiation 
field, and r denotes the thermal relaxation time of 
the electron gas. As (1.1) is linear, we may treat 
the individual Fourier components of the radiation 
field separately. The co-component of (1.1) is 



lioe u 



which has the solution 



Kl co ' 



KtI u 



are 



icpT 



€ w = 



1+W V(1 + TV) 



(1.3) 



^=tan- 1 (wr). (1.4) 



The angle <p represents the phase difference between 
the co-component of the radiation field and the 
co-component of the electron temperature. In the 
D and lower E regions, t~ 100 sec, so that the attenu- 
ation factor (l + r 2 co 2 ) _1/2 is very small except at 
audio frequencies. 

The assumption of linearity also implies that the 
variable part of the electron-molecule collision 
frequency v oc e. According to Lorentz's theory, 
k oc v. Hence, setting 



K = K + 



/- 



■/ 



we have 



K.uOCVu<X€ 



po+ v„e iat d<a } (1.5) 



(1.6) 



Finally, if the region where appreciable modulation 
transfer occurs is nearly homogeneous, we have 



XuOCKc 



(1.7) 



where X w describes the depth and phase of the trans- 
ferred modulation. From (1.4), (1.6), and (1.7), we 
have 

X^I.T-e-**; T„=(1 + t 2 u> 2 )- 1/2 , <£ w = tan" 1 (cor). 

(1.8) 

The predicted linear dependence of X u on I a , the 
predicted variation of the phase lag p^ with co, and 
the predicted variation of the coefficient T^ with co 
can all be tested separately. 

(i) Linear dependence of X^ on I u . If the disturb- 
ing wave is modulated to a depth M at a single 
frequency co and if P denotes the total power of the 
disturbing radiation, then 

7ocP[(l + iM 2 )+2Mcos co*+£M 2 cos2co*]. (1.9) 

Ratcliffe and Shaw [1958] verified that X u oc P over 
a wide range of P. Huxley et al. [1947, 1948] 
verified the predicted linear dependence of P tf on M 



2 Equation (1.1) has the same form as that governing the current in an KL 
circuit; t represents the time constant of the circuit and I the applied emf. 



and showed that the predicted quadratic dependence 
of X 2 <a on M was qualitatively in agreement with 
experiment. These experimental results show that 
a linearized description of ionospheric cross modulation 
is valid over a wide range of experimental conditions. 

(ii) The phase lag <p(o)). Ratcliffe and Shaw [1948] 
and subsequent workers found that they could 
secure good agreement between predicted and meas- 
ured phase lags over a wide range of co by choosing 
the parameter r appropriately. The required values 
of r seem to be consistent with values derived by 
extrapolating laboratory measurements [Huxley, 
1959]. 

(iii) The coefficient T a . Formula (1.8) for T« 
agrees qualitatively, but not quantitatively, with ex- 
periment [Ratcliffe and Shaw, 1948; Huxley, 1950]. 
The discrepancies between theory and experiment, 
though not large, appear to be significant. Huxley 
[1950] suggested that the two spatially separated 
regions in which cross modulation occurs in a typical 
experiment — one on the ascending, the other on the 
descending branch of the wanted ray — may have 
sufficiently different properties to invalidate the 
assumption of homogeneity underlying (1.6). Cal- 
culations based on a two-center model [Huxley, 
1950] do not, however, significantly reduce the 
discrepancies. 

The derivation of formula (1.8) actually rests not 
only on the assumptions of linearity and homogene- 
ity but also on a third major assumption: that, at 
any given point in the ionosphere, the quantities 
k, v, and t depend only on the electron temperature 
and not explicitly on the modulation frequency. 
Now, the relation between collision frequency (say) 
and electron temperature depends in general on the 
form of the velocity distribution of the electrons. If 
the distribution is maxwellian, or has any other fixed 
form depending on a single parameter, it is com- 
pletely specified by T e , so that v becomes a function 
of T e . But if the electron temperature does not 
serve to specify the distribution completely, the 
relation between v and T e may be many valued. 
In the present problem, the variable part of the 
distribution changes appreciably in a time of order 
of" 1 . Since this is short compared with the thermal 
relaxation time r there is no reason to suppose that 
the distribution remains accurately mawellian. It 
is true that if the intensity of the incident radiation 
field is sufficiently small, the velocity distribution 
of the electrons will be approximately maxwellian. 
But the phenomenon of cross modulation depends 
entirely on the variable component of the velocity 
distribution, and, as we shall see, this component 
remains nonmaxwellian even in the limit of vanishing 
field intensity. 

This paper presents a theory of cross modulation 
in which the electronic velocity-distribution func- 
tion figures explicitly. Although the present theory 
is more realistic than the macroscopic theory sketched 
above, it is oversimplified in one important respect: 
It treats collisions between electrons and molecules 
as if they were perfectly elastic. In reality, cooling 
of the electron gas in the D and E regions results 
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chiefly from the collision al excitation of molecular 
nitrogen. Basing themselves on the work of Kov- 
rizhnikh [1960], Caldirola and DeBarbieri [I960, 
L964] have recently extended the present theory to 
allow, at least approximately, for inelastic collisions. 
For the sake of simplicity, the following discus- 
sion ignores the effects of the earth's magnetic 
field. This is permissible only if all the carrier 
frequencies that figure in the discussion are much 
greater than the gyrofrequency. The modifica- 
tions required when this condition is not met are 
straightforward; they are described in the paper 
by Caldirola and DeBarbieri mentioned in the last 
paragraph. 

2. Reduction of the Boltzmann Equation 

The Boltzmann equation for an electron gas in an 
electromagnetic field has the form 



df djf e_ 

c)/ + q *dx~ i ~m 



e+^xb) 



dq 



8t 



(2.1) 



The following discussion of this equation is a gen- 
eralization of the treatment of a weakly ionized gas 
in a constant electric field given by Chapman and 
Cowling [1960, pp. 346-352]. 

We assume at first that the field is that of a 
linearly polarized plane wave propagating in the 
£ -direction: 3 



E x =Ee i{kz - pt \ 



B y =—Ee Hkz - pt) - (2.2) 
V 



The distribution function/ depends on the six vari- 
ables z, q, t. Since \ck/p\^l } the magnetic force is 
of order q/c compared with the electric force. Simi- 
larly, the second term on the left side of (2.1) is of 
order q/c compared with the first term. Neglecting 
terms of this order, we obtain in place of (2.1) 



where 



dr x du St 



F=-E 

m 



(2.3) 

(2.4) 



In this approximation, / depends on only four 
arguments: 

f=j(z;u,q;t). (2.5) 

We may expand / in Legendre polynomials of the 
argument u/q. We shall need only the first two 
terms in the expansion: 



f=f i0) (<I, t) + iFe iikz - pt) vf l} (q, t). 



(2.6) 



If the number density n of electrons is sufficiently 
small compared with the number density of molecules, 
we may neglect electron- electron and electron-ion 
collisions. The collision term in the Boltzmann 
equation, represented by 5//5£, is then a linear 
functional of /: 



fh<y,. 



Lorentz showed that 



A[«0(2)]= -v(q)ug(q)+0 



A I 



(2.7) 



(2.S) 



where m/M is the electron-molecule mass ratio, v is 
the velocity-dependent collision frequency, and g is 
an arbitrary function. Davydov [1935] derived the 
important formula 



A[0(2)]= 



mid 



= 55 !><«>( 



H 



+'[(5)1 (2 - 0) 



By inserting (2.6) into (2.3) and using the formulae 
of Lorentz and Davydov to evaluate the right side 
of the resulting equation, we obtain differential 
equations for the functions / (0) and / (1) . 4 



The equation for / (1) is 



d/(0) 



yu)_. 



1 1 df 



(0) 



(2.10) 
(2.11) 



p-\-iv q dq 
The left side of the equation for / (0) has the form 

(2.12) 



d/ (0) / d \ 



Here the brackets indicate a double average: over 
direction in velocity space (to eliminate terms con- 
taining Legendre polynomials P h (u/q) with F>0), 
and over time (to eliminate radio-frequency fluctua- 
tions in y (0) , which are of no physical interest). 
Since we need to preserve audio-frequency variations 
of / (0) , the time-averaging must be over an interval 
that is short compared with the reciprocal of the 
highest modulation frequency of interest as well as 
long compared with the reciprocal of the carrier 
frequency p; thus we must have /;»co. Applica- 
tion of the well-known rule 



(Re Ae ipt • Re Be ivt ) t =\ Re AB, (2.13) 



3 We adopt the convention that when a real quantity, such as E x or/, is repre- 
sented by a complex expression, as in (2.2) and (2.6), the real part of this expression 
is to be understood. Note, however, that the quantities E, F, fc,/(°> ,etc, are all 
complex, and that both the real and imaginary parts of an equation like (2.11) 
are significant. 



4 The differential equation for/(*>, the coefficient of Pk(u/q) in the Legendre 
expansion of/, involves only the fund ions/' " with j<k, so that the sequence of 
approximations may be terminated al any point. By contrast, the equations 
for the moments of/ of order k involve moments of order A;+l. 
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where A denotes the complex conjugate of A, gives given by 

1*/ \F,f tT, 

X instead of by (2.17). 

=-Re/?f ^ V 5/^d/ (0) \\ 
2 \ \jP+W S^vz &2 / J 3. The Dispersion Relation 

1 j/ . j~ l d/ (0) , w 2 5 /l 5/ (0) \"l (cy*iA\ ^ ne e l e °t r ic current density produced by the 

_ 2^H7 2 ' ' U~d<r 7dg\2^F/J radiation field (2.2) is given by 

(F&tiF«n) -<;'-JMt*&} (3A) 

__1 v IFI'P^Vg a ^ ld -H1 J,=J.=0. (3.1') 

2p 2 +x 2 ' ' U dg^3 a d2\g £><? /J 

From Maxwell's equation, 



Since 



averaging over direction in velocity space gives 



<Sp 2 + v 2] ] q 2 dqV dq J curl B=- E+4ttJ> (3.2) 

(2.15) 

and (3.1) we obtain 

Combining this result with (2.9) we have, finally, , „ . , 

(ck) 2 = p 2_^P Cfm d (l*_£_\ 

df w m 1 d / 3 ( Sw ,eaf»l\ n /oi^ mJ V3p+i"/ 

tr- 2?5i {-<4-' , " + 7^J>^ < 216 > M . f ,.^ 



where 

~ 1M IFP 



We may immediately generalize these results to an J n order *° wr j te thi f relation in a more compact 
arbitrary radiation field. The electric vector at a J?™ 1 ™e introduce the plasma frequency p , de- 
given point may be written in the form fined by 



_ 4wne 2 



E=SE„ (2 . 18) fi=~- (3.4) 



i i -Tj, , ^ , ,i and use the abbreviation 
where each Courier component % represents the 

electric field of a linearlv polarized plane wave. In 1 r T4 "1 

place of (2.6) we have <$)=± \f 0) d\ ~q 3 ^(q) \> (3.5) 

J—J° (j^+Ij^w'V/j (<Z>0 ex P (^Ky-x— Pjt) where n denotes the electron density. The disper- 

ro ^Q^ s ^ on re l a ^i° n (3.3) then takes the form 

and in place of (2.10), {c kf=f-fJj^)+i P l(j^-\ (3.6) 



/)" = X- 1 ^— (2-20) 

Jpj\ %v( l 0C L The absorption coefficient k is defined by 

Finally, the function / (0) satisfies (2.16) with 6 Jc=Jc +%iK, (3.7) 
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where k and k are both real. The following approxi- 
mate formulae, valid under the conditions stated, 
are often convenient: 



ck= iA/pl\ ( Po<<p ), 

VxP\~ v / 



(3.8) 



CK = 



(f) 2 « (Po«P,»«P). (3.9) 



The results of this section were first obtained by 
Lorentz. 

4. Absorption of Unmodulated Radio Waves 

When co = 0, (2.16) reduces to an ordinary first- 
order differential equation, whose solution is 



f Q) =7iCexp\ 



-r 



e-w» 



(4.1) 



where 9 is given by (2.21) and the constant C is 
determined by the normalization condition 



4 "/> 



) q 2 dq=n. 



(4.2) 



Formula (4.1) is valid for all values of the carrier 
frequencies pj occuring in formula (2.21) for 9. 
We recall that / (0) represents the isotropic part of 
the distribution function averaged over a time inter- 
val long compared with all the characteristic times 
pj l . If the radiation field has only a single Fourier 
component, / (0) represents the isotropic part of the 
distribution function averaged over a single period. 
Formula (4.1) remains valid in the limit p = 
(E= const.). In this case 



9 = 



1M\F\ 2 kT m 
6 m v 2 (q) m 



(p=0) 



(4.3) 



and (4.1) coincides with a formula derived by Chap- 
man and Cowling [1960; eq (13), p. 350]. 

When the carrier frequency is very large, 9 is 
given by 



9= 



\M\F\ 2 .kT r 



6 m p 2 



+- 



m 



(P»v). 



/(°) is accordingly maxwellian, the electron tempera- 
ture being given by Lorentz 's formula, 



kl e — kl m -\- 



\ M\F\; 
6 p 2 



(4.5) 



If the molecules are maxwellian (collision frequency 
independent of velocity), the distribution function 
is maxwellian for all values of the carrier frequency. 
Otherwise the distribution function departs markedly 



If the molecules are rigid elastic spheres (v = q/l, 
1= const.) the integral in (4.1) may be evaluated 
explicitly; one obtains 



fO) = C f (q 2 +q 2 +A) A e-^, 



where 



mq l 
~'2kf m 



- 2 _ml 2 p 2 

P ~2kT7n 



A= 



M\F\ 2 



2kT m 



(4.6) 

(4.7) 



■(0) 



5. The Variable Part of / 



In this section we derive the equations governing 
small departures from the steady-state velocity 
distributions discussed in the preceding section. 

Let 



e(<Zi*)=3o(<z)[H-%,*)]. 



(5.1) 



The function 9 coincides with the time-independent 
function that was called 9 in the preceding section. 
We define new dimensionless variables x, X, <p: 



dx- 



_gdz. x--^ 2 - 



9 



2 9 



^(»,*;=4im- i r f 2y w, (2 l *)d2. (5.2) 

Jo 



Note that 



dx_dX d9 () 
X~~X^^ 



(5.3) 



The function <p is the probability distribution function 
of the dimensionless velocity variable x; 



ip (a,t)=Vr{x<a}. 



It follows that 



*((),*) =o, *>0»,*) = i. 



(5.4) 



(5.5) 



The partial derivative d<p/dx is the probability 
density associated with x. In terms of the new 
variables (2.16) takes the form 



dcp 2m 



(4.4) dt M 



l\d* 
2/ dx 






So far we have not made any assumptions about 
the magnitude of the function 0(q,t). In order to 
reduce (5.6) to an ordinary differential equation we 
now set 



<p=<Po(x)+<Pi(x,t) 



(5.7) 



and assume that 6 and <p { are so small that their 



from the maxwellian form when 9>> kT m /m and product may be neglected in (5.6). This will be 
p~v. true if at least one of the conditions 9«k2 7 , r /m, 
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744-758—65- 



M<<1, where M denotes the modulation depth, 

is satisfied. Having linearized (5.6), we can deal 

separately with the Fourier components of the s0 that (5.12 )becomes 

radiation field, just as in the elementary theory of 

section 1. Let 



X—x, o =const., 6 u =const., (5.16) 



S(x i t)=B u (x)e i 



*>iCM) =*>«(*)*'"'• (5-8) 



It is convenient to replace the modulation frequency 
co by a dimensionless variable and to separate out 
the velocity dependence of the function v(x). We 
accordingly define 



V(x) = 



v(x) 



v= I v(x)d<p 



2(m/M)v 
Finally, since ^ satisfies the equation 



(5.9) 



(5.10) 



™>u:m x -2)is^ t«» 



we have 



dx 



=ce 3 /2 x 1 ' 2 e - x , 



(5.12) 



where C is determined by the normalization con- 
dition <p ( oo ) = 1 . Hence 



-X^=-C(Q Xy< 2 e- 



(5.13) 



Omitting second-order terms in (5.6) and using 
(5.9), (5.10), and (5.13), we obtain 



'dx 



( e " i- &7 + (^ 2jdx 

+ iar J( p u =d 0) Ce 3 /2 X 3/2 e- x . (5.14) 



Since <p and <p both satisfy the boundary conditions 
(5.5), <p w must satisfy the boundary conditions 



^(0)=^(oo)=0. 



(5.15) 



Equations (5.14) and (5.15) together with the 
dispersion relation (3.6) represent the formal solution 
of our problem. The functions O and 9 are defined 
by the radiation field. Given 6 , one finds the func- 
tion x(q) and X(q) from (5.2). Equation (5.14) 
with the boundary conditions (5.15) may be inte- 
grated numerically by the method described in the 
appendix. Finally, knowing <p u , one can calculate 
the variable part of the absorption coefficient from 
the dispersion relation. 

Equation (5.14) assumes a simpler form when 
G~kT m /m, so that the unperturbed distribution 
function is nearly maxwellian. We may then write 



Setting 



<^o=JL x 1/2 e- x 
dx ^ T 

2 

<Pa,(x) = --=d (a y(x), 

■yjw 



we obtain in place of (5.14) 

xy" + (x—%)y'—iar)y= 

The boundary conditions are 

2/(0)=2/(°°)=0. 



(5.17) 
(5.18) 

(5.19) 
(5.20) 



In the numerical work described m the next section 
the function rj(x) is taken to have the form 



Vtt 

6. Some Numerical Results 



(5.21) 



Miss Cara Joy Hughes employed the IBM 7094 
computer at the Harvard University Computing 
Center to integrate (5.19) with the boundary con- 
ditions (5.20) and with rj(x) given by (5.21). In- 
tegrations were carried out for every pair of param- 
eter values (a,r) with a in the set (0.1, 0.2, 0.5, 1, 2, 
5, 10, 15) and r in the set (-.5, .5, 1, 1.5). The 
integration procedure is described in the appendix. 

For r=0 the solution of (5.19) is 



y= 



1 + ia 



(r=0). 



(6.1) 



The complete distribution function f (0) corresponding 
to this solution is maxwellian with a variable electron 
temperature. The factor (1+ia) -1 in (6.1) cor- 
responds to the factor (l+iwr) -1 in (1.4). 

Figures 1 and 2 show the modulus and argument of 
the function 



y=(l+ia)y 



(6.2) 



for r= —0.5, 1.5 and for a=2, 5. In both figures the 
curves corresponding to r=— 0.5 and r=1.5 differ 
markedly from the curve for the maxwellian case 
r=0. For r=— 0.5 the peak of the function \y(x)\ 
occurs at a smaller value of x than for r=0, while for 
r=1.5 it occurs at a larger value. In the four cases 
with r^0 the function arg y changes sign near 
x=\.b. For large positive values of x, arg y is 
positive for r=1.5, negative for r=— 0.5. 

According to (3.9), the absorption coefficient k 
varies directly as the mean collision frequency 
(the mean being defined by (3.5)) if the carrier 
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Figure 1. The function |y| (see eqs 6.2 and 5.18) for 
representative values of r and a. 



frequency is sufficiently high. By (3.5), (5.2), 
(5.7), (5.8), (5.18), and (5.21), 

(p a )oz f (I ^ if 2 d(<fv)ozC2r+Z) fd<p„x' 

f* co /• 00 

oc - (2r+3)r ^ w x r ~V/xoc - (2r+3)rda yx r ~ l dx. 

Jo Jo 

(6.3) 

Just as in section 1 , the assumptions of linearity and 
homogeneity imply that X u oc k u (see (1.7)). Hence 



X„oc(2r+3)rd ( 



yx r ^x. 
Jo 



(0.4) 




Note that Xu vanishes for maxwellian molecules 
(r=0). This means that the elementary theory 
becomes rigorously valid in the limit when the 
phenomenon it describes disappears. 



By analogy with (1.8), we write 

/» GO 

yx r - l dx=T(a)e- tria) - 



(6.5) 



Figure 2. The function arg y (see eqs 6.2 and 5.18) for 
representative values of r and a. 



The functions </>(o:) and T(a) are shown in figures 3 
and 4 for r=-0.5, 0, 0.5, 1, 1.5. For r>0 the 
curves </>(a) have roughly the same shape, (im- 
pressing the horizontal scale of the curve for r = 0.5 
by 25 percent would bring it into near coincidence 
with the curve for r=0, but for r=l, 1.5, the scaling 
factor increases markedly with x. For r=— 0.5, the 
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a 

Figure 3. Phase lag of the transferred modulation (see eq 6.5) as a function of modulation frequency for various 
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collision laws. 
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Figure 4. Amplitude of the transferred modulation (see eq 6.5) as a function of modulation frequency for various collision laws. 



phase lag at first increases more rapidly with 
increasing modulation frequency than in the max- 
wellian case, then more slowly. The shapes of the 
curves for r=— 0.5 and r=0 are entirely different. 

If one were to use the phase-lag curve for r=0 to 
analyze experimental results relating to a hypo- 
thetical gas in which the electron-molecule collisions 
were elastic and were characterized by a velocity- 
independent free path (the case r=0.5), one would 
over estimate the mean collision frequency v by 25 
percent. If a higher value of r were appropriate, 
the error could be much greater. 

Turning now to the coefficient of transferred 
modulation T{a) (fig. 4), we see that the curvature 
of the function T(a) increases with increasing r over 
the entire range — 0. 5<r< 1.5. Compressing the 
horizontal scale of the curve for r= 0.5 by 25 percent 
would make it fall off more steeply at small and 



moderate values of a than the curve for r = 0. In 
general, one cannot devise a horizontal scale trans- 
formation that, for a given value of r^O, will make 
both T(a) and <p(a) assume the forms appropriate 
to the case r=0. 

To sum up, the numerical calculations show that 
significant departures from the predictions of the 
elementary theory may be expected. Moreover, 
values of the mean collision frequency derived by 
using the elementary theory to analyze cross- 
modulation data may be significantly in error. On 
the positive side, the sensitivity of the predictions 
to changes in the collision law suggests that further 
theoretical and experimental refinements could 
ultimately lead to an accurate experimental determi- 
nation not only of the collision frequency but of the 
form of the collision law. 
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7. Appendix 

Cam Joy Hughes 

Division of Engineering and Applied Physics, Harvard University 

The following paragraphs describe the method used to integrate (5.19). The same 
method can be applied to the more general equation (5.14). 



Consider the inhomogeneous second-order linear where L and R have the forms 
differential equation 



a{x)y' r + b{x)y' + c{x)y=d{x). 



(A.l) 



We may approximate (A.l) by the set of coupled 
difference equations 

aA 2 lh'rbAy i + e i y i = d i (i=l, . . ., n-\) (A.2) 
where 



a i = a(x i ), etc. 
Xi=x +ih 



(A.3) 

(A.4) 

(A.5) 



1 
l 2 1 
I, 1 



• ln- 



R= 



(A.2) is a set of (/? — 1) equations For the (ft+1) 
variables y ,yi, ■ . ., y n > r rtie two boundary condi- 
tions provide two additional equations, so that in 
general the set of difference equations together with 
the boundary conditions have a unique solution. 
In the problem at hand, y and ?/„ arc 1 given: 
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n o . . 


. 





k 2 r 2 . . 


. 





o h n • 


. 







■ K 



(A. 10) 



2/o=2/n=0. 

The equations (A.2) thus have the form 



The coefficients /,, /•,, r x must satisfy the equations 

(A.6) 7rill== ] Ch mi 2 =n, m <i< _i=Z|ft i _i, ^«,i=^^-i+^, 

m M+1 =r, (2<i<n-l) (A.ll) 



/ 2a,,_, \ ,/a„_2 &«-2\ , 



In matrix notation, 



MY=D 



(A.7) 



(A.S) 



where M is a tridiagonal square matrix of order 
(n— 1), and F, Z? are column matrices. 
To solve (A. 8), we write Min the form 



which may clearly be solved in serial order for k i} r h 
l 2 , k 2 , r 2 , etc. Thus L and R are completely specified 
byM. 
Let 

RY=U. (A. 12) 

Then (A.S) becomes 

7X r =L> (A.13) 

or, again equating matrix elements, 

Wi=di, J«u,-i+tt«=eJ, (2</<//-l) (A.13') 

which can be solved serially for the u t . Having 
found U, we then determine Y from (A. 12), giving 

k i y i +r i y i+1 =u i (l<i<n—2) 

kn-iy f i-l = 'Un-l (A.12') 



M=LR 



which can be solved in reverse serial order, beginning 
with the last equation in the set. The solution is 
(A. 9) now complete. 
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One can estimate the accuracy of the procedure de- 
scribed, as applied to the differential equation (5.19), 
by comparing numerical solutions of this equation for 
tj = 1 with the exact solutions (6.1). The numerical 
solutions were carried out for a grid spacing h =0.05 
and with the approximate boundary condition 
y(x=25) =0 in place of y(x= «) =0. 
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